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Abstract. Using complex analysis techniques we obtain precise asymptotic 
approximations for the kernels corresponding to the symmetric a-stable pro- 
cesses and their fractional derivatives. We apply our method to general Levy 
processes whose characteristic functions are radial and satisfy some regularity 
and size conditions. 



1. Introduction 

A stochastic process {Xt : t > 0} defined on a probability space (17, J", P) is said 
to be a Levy process if Xq = almost surely, for any s,t > 0, the increment 
Xt+s — Xt is independent of X^ for < u < t and has the same law as Xg, and Xt 
is almost surely right continuous with left limits. By the definition, its characteristic 
function is given by 

E(exp(ze-XO) -exp(-i7y(0), 
where 7y is a Levy symbol. The symbol i] plays a central role. Actually, it is possible 
to characterize all Levy processes by their Levy symbol rj, in the sense that, one can 
construct a Levy process from any symbol of the form given in the Levy-Khinchine 
representation 

viO = ^ - • e - / (e'^-" - 1 - ^^ xl{|,|<i} {x)Hdx), 

where A is a nonnegative-definite symmetric matrix, b G M.'^, and the Levy measure 
v is a positive Borel measure satisfying 

/ A < oo. 

jR<i\{0} 

Here, we use the notation a Ab = min {a, b}. We note that the index of rj given by 

a = inf|A>0: lim ^ = ol 

should satisfy < a < 2 (see Blumenthal-Getoor [BCiGl]). 

Let us mention some of the fundamental symbols which are radial. Brownian 
motion is a basic example which lies at the source of the theory of stochastic pro- 
cesses. A symmetric a-stable process is another important example whose symbol 
is given by 

77(0 = |?r, 0<a<2. 
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We note that the infinitesimal generator of a symmetric a-stable process Xt is the 
fractional Laplacian — (— A)"/^ that can be written in the following principal value 
integral 

(-Ar/M-) = Cp.v. / "P^^dy 
^ ' ^ ' V \x-y\''+" 

for some constant C . Recently there has been many interests from the theory of 
probability, PDE and SPDE to study such fractional Laplacians (see e.g., Caffarelli- 
Silvestre [CS07], Chen-Kim-Song [CKSIO], Bogdan-Sztonyk [BS07]). We now men- 
tion one more example. A relativistic stable processes whose symbol is given by 

77(0 = (|ei' + m2)"/2-m" 

is an example of radial symbol, but it has no scaling property for m 7^ 0. Note 
that when m = 0, this coincides to the symmetric a-stable process. Recently, 
sharp kernel estimates for relativistic stable processes in open sets are obtained by 
Chen-Kim-Song [CKS12]. 

Analytic tools such as semigroup, resolvent and infinitesimal generator have deep 
connections with the probabilistic properties such as transience and recurrence, 
asymptotic behavior at infinity, etc. The family of convolution operators indexed 
by t > and given by 

Ptf{x) = f f{x + y)¥{Xt e dy) 
is a Markov semigroup. Taking the Fourier transform, we write 

^{Ptfrn = E / e^^--f{Xt + x)dx = 
Thus, it is very important to obtain sharp estimates for kernels Pt{x) given by 
Pt{x) ^ {2tt)-'^ I e~'=^-«e-*"(-«)df, /3 > 0. 



For more information about Levy processes, see Sato's monograph [Sat99]. 

In this article, we propose a method to obtain the asymptotic expansion for 
< a < 2 

P^{x) = {2tt)-'^ ( e-"-«e-*l«l°dC 
and its fractional derivatives for /3 > 



Interestingly, for the latter case the decay order does not vary continuously on 
/? S K. We also apply our method to the general radial symbol 77 which has some 
regularity conditions similar to the standard symbol class in the theory of pseudo- 
differential operators. There are many results about asymptotic behavior of stable 
laws. The asymptotic expansion of one-dimensional stable laws was obtained by 
Linnik [Lin-54], Skorohod [Skc).")4], Zolotarev [Z0I8G], and others. There are few 
result for the asymptotic expansion of stable laws in dimension d greater than 1. 
We should mention the work of Kolokoltsov [KolOO] that is similar to ours. He used 
special functions heavily such as Whittaker functions and then he generalized from 
his results to a different direction compared to ours. We try to use effectively the 
deep connection between the decay of the kernel and the singularities of the Mellin 
transform of the symbol. Our method of proof is elementary and direct. It is based 
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on the very basic but decisive complex analysis technique, that is Cauchy's Residue 
Theorem with shifting contour integration. As a byproduct, c^xact constants are 
also obtained. The Bessel functions naturally arise from the spherical symmetry 
and we will use some of its basic properties. 

The organization of this paper is as follows. In Section 2 we introduce the 
Mellin transform and its inversion formula and give an important example which 
is the gamma function. Some of the basic properties of the gamma functions and 
Bessel functions are also presented for our purpose. In Section 3 we prove complete 
aymptotics of the kernel corresponding to the symmetric a-stable processes via 
Mellin's trasform and give some corollaries. The proofs in this section show that 
the asymptotic behavior of an oscillatory integral is closely connected with the 
singularities of the Mellin transform of that integral. In Section 4, by refining 
the same techniques already developed in the previous section, we investigate the 
asymptotic behavior of the kernels corresponding to more general radial symbol r] 
under some weak conditions. 

2. Preliminaries 

Throughout this article we shall be working on R*^ with d > 2. Given an inte- 
grable function f{x) on R"^, we define its Fourier transform Tf{^) by 

-^/(O = / e'^--f{x)dx, 

where ^ G R'' and ^ ■ x = Y^'jLi^j^j- Similarly the inverse Fourier transform is 
defined by 

Moreover, if / and T f are uniformly continuous and integrable on R'', then the 
Fourier inversion formula holds, that is, 

f{x)=F-\Tf){x). 

More generally, the space of integrable functions on is contained in the space of 
finite complex measures on with the total variation norm via the identification 
fdx = dji. Thus we can generalize the definition of Fourier transform to measures 
via 

= I e*«-d/x(x). 

Given a function g{r) on [0, oo), we define its Mellin transform M.g{z) by 

/•OO 

Mg{z) = / g{r)r''~^dr. 
Jo 

This integral defines a holomorphic function on a strip in the complex plane. It 
is closely connected to the Fourier transform. Actually, it can be expressed as a 
Fourier transform. To see this, let 2: = a + i6 and then change variables r = to 
obtain 

/oo 
g{e-^)e-''='e-''"'dx. 
-00 
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Mellin's inversion formula also follows from the Fourier inversion formula. To see 
this, if 5((e-^)e-"^ e L^{R,dx) and Mg{a + ib) e L^{R,db), then by the Fourier 
inversion formula 

1 f°° 

5(e-^)e-°^ = — / e-'''''Mg{a - ib)db. 
Hence we get the Mellin inversion formula 
(1) g{r) = — / Mg{z)r-'dz 

by the change of variables. From now, we shall use the notation 

f{z)dz = / fiz)dz. 



' (c) «/ c—ioo 

A basic example is the gamma function which is defined by the Mellin transform 
of e~*. It is initially defined for > by the absolutely convergent integral 

f-OC 

r{z)= / e-H'-^dt. 



We now present some of basic properties of the gamma function. Integrating by 
parts gives the functional equation r(2: + 1) = zr{z) and this relation produces a 
meromorphic continuation. We notice that there is another method to obtain the 
meromorphic continuation, splitting the integral and expanding e~* in a power 
series also produces 



r{z)= y^-^t'-'dt+ e-H'-'dt 
= E + / e-H^-'dt. 

The series defines a meromorphic function with simple poles at z = 0,-1,-2,... 
and the integral on the right defines an entire function. Thus, the gamma function 
has an analytic continuation to a meromorphic function on C whose singularities 
are simple poles and the residue at z = —n is (— l)"/n!. From Euler's reflection 
formula 

r(z)r(i-z) = ^, 

smTTz 

l/r(z) is an entire function with simple zeros at z = 0,-1, —2, . . . and it vanishes 
nowhere else. 

The Bessel function of the first kind of order v can be defined by a numbcir of 
difii'erent ways. One of the usual ways is the Poisson representation formula which 
is given by 

r(j. + i/2)r(i/2)y_i'^ 

where ^{v) > —1/2 and r > 0. Finally, we note that the Mellin transform of 
Jv{r) is given by 



(2) I 
Jo 



r-'J^ry-'dr 



2^—'TC^z) 
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for V > -1/2 and for < 5R(z) < + 3/2. We shall use this identity frequently. 
For more information about Bessel functions, see Watson's monograph [Wat95]. 

3. Symmetric q;-stable processes 
We begin by considering the classical heat kernel 



We note that it has the following scaling property 

This reflects selfsimilarity of the Brownian motion, that is, any change of time scale 
for the Brownian motion has the same effect as some change of spatial scale. Thus 
we may assume t — 1 and use Fubini's theorem to obtain 

/>oo 

Pi(x) = (271)-'^ H / e-^^"«"e-«'d^„. 

2 . 

The function is holomorphic in the entire complex plane and has rapid decay 
as — > oo. Hence by Cauchy's theorem we can shift the contour integration 
to yield 

-OO J — oo 

Therefore 

P,(x) = (47rt)-'*/2e-l^l'/4*. 
We now consider the symmetric a-stable processes 

-d f -tier 



(x) = (271)-" / e-"-^e-'l^l dC. 

The same argument breaks down in this case because its symbol has no multiplica- 
tive structure and has singularity at the origin. However, to prove the following 
theorem, we shall effectively use the above fundamental technique, that is, to trans- 
form the multi-dimensional integral to the contour integration and then to evaluate 
the integral by shifting the contour. 

Theorem 1. Let d > 2 and < a < 2. Then P^{x) has an integral representation 
for (d — l)/2 < c < d and an asymptotic expansion 

(4) ^ Ed - ^^(W2)) \y \ \l {t~^'-\x\r^—. 

n—l ^ 2 / 

Proof. We may assume that t — 1 because the kernel has selfsimilarity 

(5) (x) = t-'^/"P^{t~^/°'x). 
Claim 1: 

/■oo 

(6) P^{x) = (2^)-'^/2|x|-''/2+i / J,/^_,{\x\r)e-^''r''/^dr. 
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Proof of Claim 1: In polar coordinates, 



(27r)-'' / J-fj 
Jq 



where da is surface measure on S''^^. Now (6) is a consequence of the well-known 
identity 

J-a(x) = (27r)'^/2|xr''/2+ij,/2_i(|x|), 

where J^, is the Bessel function of the first kind of order v. For the proof and 
much more general information about Fourier transforms of measures supported on 
surfaces, we refer to Stein's monograph [ ki , ,]. 
Claim 2: For (d - l)/2 < c < d 

Proof of Claim 2: We recall the gamma function 

r{z) = / e-W-^dr, 



which is the Mellin transform of e ' . From Mellin's inversion formula (1) we have 
for c > 

1 



e-'' ^— / T(z)r-''dz. 
Thus we obtain, by rescaling, the key identity 
(8) e"'" = — / a-^r{a-^z)r-'dz. 

Inserting (8) to (6) we see that making the change of variables yields 

P^{x)^{2tt)-'^/^\x\-^'^+^ r Jd/2-i{\x\r)^. f a-^T{a-^z)r-'dzr'^'^dr 



= (27r)-'^/2|^|-d/2+i 1 r ^-ir(c,-iz) r J^/^_,{\x\r)r-^+'"^drdz 

= {2t:)""M-^^. [ a-'Tia-'z) H J^/^-^iry-^+^'^dM^dz. 

The identity (2) gives for (d - l)/2 < c < d 

'Jd/2-iKr}r dr - ^^^^^ . 

Therefore the integral representation (7) follows. 
Claim 3: 



oo 



(9) = E(i - i^(W2)) ^j \\L 

n— 1 ^ 2 / 
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Proof of Claim 3: To obtain this series representation we consider first the integrand 
in (7) 

(10) ni) 

We notice that r(^) is a meromorphic function with simple poles at z = 0, —a, —2a, . 
and that l/r(|) is an entire function with simple zeros at z = 0, —2, —4, .... So if 
na/2 ^ Z, then z — —na is the pole of the integrand (10). 

We define the contour 7^ consists of a rectangle with vertices c~iR,c + iR, — c + 
iR, —c — iR with the counterclockwise orientation. Now, we recall the following 
version of Stirling's formula for the gamma function. If a < m < & and \v\ — >■ 00, 
then 

(11) \T{u + i«)| = V2^\vr^'^ ( ~ 1 1^1) (1 + , 

where the constant implied by O depends only on a and b. By Stirling's formula 
the integrand (10) has rapid decay as |9(z)| — 00. Hence the two integrals over the 
horizontal segment on top and bottom tend to as i? tends to infinity. By Cauchy's 
residue theorem we can shift the vertical line integration to the left c — > — c and 
the remaining integral tends to zero as \x\ — > 00. Therefore by letting c -> —00 we 
obtain the desired expansion (9). 

Finally, using the relation (5) we obtain (3) and (4) from (7) and (9). □ 

Remark 1 . This proof shows that the asymptotic behavior of an oscillatory integral 
is closely connected with the singularities of the Mellin transform of that integral. 

Remark 2. For the case a — 1 the asymptotic expansion (4) coincides with Poisson's 
kernel. In fact, we have 

^<i/2|^|d+i A. (2n + l)!r(-n-i) ' ' ' 

A calculation gives 

(2n + l)!r(-n - 1/2) = 2^"+^nl{-l)''+^r{l/2) 

and so we have 



-n( d+1 \ 00 , d+1 ' 



r(^) 



_^^(i+Nr^- 

TT 2 

Remark 3. For the exceptional case a — 2 the integral identity (3) is still true. So 
we have 
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By Cauchy's residue theorem we can shift the contour integration to the right to 
obtain 

' ' n=0 

= (47r)-'^/2e-l"l'/4. 

We denote f < g ii f < eg ior some constant c > and denote f ~ g it f ^ g 
and g f ■ The foUowing corollary is a direct consequence. 

Corollary 1. Let d>2 and < a < 2. 
Proof. For > 1 we have from (4) 

The result follows from (5). □ 
We now consider the symbol 

with < a < (3 < 2 that has no scaling property. For this case we have the 
following interesting corollary which shows the decay order depends mostly on a 
for fixed time t but the scaling relation depends on time. In the next section we 
study the general symbol of the form 

under some assumptions on r]i. 

Corollary 2. Let d > 2, < a < ^ < 2 and 

K,{x) = {2n)-^ f e— « exp(-t[|Sr + \if])dt 

Then 

(Kt{x)<t-'^/^{l + t-^/^\x\)-<^'^+"^ ift<l 
[Kt{x) < ift>l. 

Proof. Note that 

Kt{x)^ I Pt{x-y)Pt{y)dy. 



First we assume t < 1 and consider the cases t ^^^\x\ < 1 and t ^/^\x\ > 1. In 
the case t~"'^/^|a;| < 1 we have 



(l + i-i//3|a;|)d+a- 
Therefore 

2d+a-^-d/fi 



Ktix) 



< 



{i + t-yf^x\y+'- 



[ P^{x-v)dy. 
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In the case t ^^^\x\ > 1 we split the mtegral to obtain 
Ktix)< [ Pr{x-y)Pfiy)dy 

J\y-x\>\x\/2 

'\y\>M/2 

^ , .-M.^^^^,^. \ Pf {y)dy 



Since t-'^/'^\x\ > 1 and t'^'°'/l^ < 1, we have 



P,^{x-y)dy. 



t^d+a-^—d/ a 2^^^t ^^d+f^-^—d/l^+l — a/lii ^d+a-^—d/ (5 

< ; = -. < 



(l + i-i/"|a;|)'^+" - {t-^/f^x\Y+» ' {1 + t-^ P\x\)d+'^ 

and the required estimate foUows. 

Now we assume t > 1 and consider the cases < 1 and > 1. In 

the case < 1 we have 

nd+a-t—d/a 

P°'(r - lA < i^'^/" < 

*^ - [i + t-y»\x\Y+^ 

Therefore 

Kdx) < 7 , / Pf(y)dy. 

In the case t~^/"|a;| > 1 we spht the integral as before to obtain 

od+aj.—d/a r 

KAx) < 7 , / Pf{y)dy 



2d+p^-d/l3 



Since > 1 and ^i-/^/" < 1, we have 



/ P^{x-y)dy. 

J\v\>\x\/2 



2d+l3^-d/p 2'^+'^t ^d+l3-^-d/a+l-l3/a ^d+P^-d/a 
< , , , ,. < ; ^ < 



{i + t-^/p\x\Y+p - ixi'^+i^ - {i + t-^/^\x\y+p - + 

and the required estimate follows. □ 

Now we are using multi-index notation here and note that formally for multi- 
index B 



More generally, the fractional derivatives of Pt{x) is defined at least formally by 



Actually, this is the Riesz potential of the probability density function P^{x). Some- 
times it is necessary to estimate the fractional derivatives of Pt{x) especially when 
one studies a certain type of stochastic partial differential equations. For example, 
in Chang-Lee [CL12], Sobolev and Besov estimates for an SPDE with a fractional 
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Lapalacian operator was obtained from the decay estimates of P"(x) and it's deriva- 
tives. The following corollary is rather surprising. 

Corollary 3. Let d > 2, < a < 2 and /3 > 0. Then {-/\Y/'^P^{x) has an 
integral representation 

j.-{d+l3)/a 1 r z \-n ( d+P-z z 

27rzi(,) r(^) ^ ' 

// /3 7^ 0, 2, 4, . . then as \x\ — > oo the first order approximation becomes 

9/3-2 -pf d+P \ 

If (5 — 0,2,4, . . ., then as \x\ — > oo the first order approximation becomes 

Proof. We may assume that t = 1 because the kernel has the following scaling 
property 

(_A)'3/2pf (x) = i-(''+'5)/"(-A)'3/2p^"(t-l/"2;). 
By the same way to derive (6) we also have 

/•oo 

i-Af^^P^ix) = (2^)-'*/2|a;r'^/2+i / j^/2-i(|2:|r)e-'^\'^/2+^dr. 

Using the identity (2) we have 

{-Af/'P^{x) = / ^"^^ Ixj-'^-'^+^d^. 

If (3 is not an even integer, then is the first pole of the integrand. Similarly, we 
shift the line of integral for si from ci to oo so that we obtain 

as |a;| — oo. 

If /3 is an even integer, then —a is the first pole of the integrand. In order to get 
an asymptotic expansion near infinity, we shift the line of integral for z from c to 
00 so that we obtain 

as |a;| — >■ oo. □ 

Remark 4. An examination of the proof shows complete asymptotic expansions can 
be obtained. 

Here is a direct consequence. 

Corollary 4. Let d > 2, < a < 2 and /3 > 0. If ^ ^ 0,2, A, ... , then 
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///3 = 0,2,4, then 



|2,|d+/3+a 

Proof. Since the function |^|^e~l^l is integrable on R'^ for all /3 > 0, the function 
{—A)^^^p{l, x) is bounded and smooth on (0,oo) x M'' by the Riemann-Lebesgue 
lemma. □ 

Remark 5. In Kim-Kim [KK12], the authors obtained an upper bound 

\{-Af/^P^{x)\ < \x\-'^'^+'^\ 

This was crucial for them to prove their generalization of a Littlwood-Paley operator 
for the fractional Laplacian. It was remarked that they thought this upper bound 
would not be sharp by considering the case /3 = 0. However, this corollary says 
that the order of decay is optimal. 

4. Radial Symbols 

In this section we refine our argument in the previous section and give a few of 
applications which are generalizations of the previous results. We shall study the 
asymptotics of the kernels corresponding to the radial symbols 77(|^|) in M'', d > 2. 
From now we consider the following general kernel defined by the integral 



Kf'{t,x) ^{2Tr)-^ e-'"-«|C|''e-*''(l«l)d^ for /3 > 0. 
By the same way to derive (6), we obtain 

noo 

(12) if'3(i,a;) = (27r)-''/2|^|-<i/2+i / j^^^_^(^\x\ry^^+^e-''^^'''>dr. 



Here is our first general result. The key point is to find a Mellin's inversion formula 
for the symbol e"*''^''^ and then to obtain a meromorphic continuation which can 
be obtained by performing an integration by parts. 

Theorem 2. Suppose that rj{r) satisfies 



(13) liminf- 

1 — >oo log r 

and for some < a < 2 and k > (d + 3)/2 + /3 

(14) sup sup |r-"+™L>'"?7(r)| ^ A < oo, 

0<m<k r 

where D denotes the derivative with respect to r. If /S ^ 0,2, A, ... , then for each 
t > 



K^{t,x) 



as \x\ — ^ oo. 



^d/2r(^) \x\d+P 

Proof. Let Mt{z) denote the Mellin transform 

/"OO 

Mt{z) = / e^*''(''V"-idr. 
This integral converges absolutely for 'R{z) > by the condition (13). 
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Claim 1: For 3?(z) > 
(15) M,(z)= ^^J^^^^ Af,^(z), 

where 

nOO 







Proof of Claim 1: If each of the boundary terms vanishes, then the repeated inte- 
gration by parts yields 











z + l)zJo ^ ' 



n^+k) Jo ^ ' 

To see that the boundary terms vanish, we use the condition (14) to obtain 
|^£)(e-*''('-))| = |tri:>77(r)|e~*"(''' < Atr^e-*''^'') 

and 

|^2^2^g-t^(r)>)| < |ij.2i:)2^(r)|e-*''('') + |iV^D77(r)|2e-*''('') < (Air"+(Air")2)e-*''(''). 
Inductively, we can obtain for 1 < m < fc 

(16) |r'"i:)"(e-*"(''))| < (Atr" + • • • + (Atr")™)e-*''(''^ 

Consequently we have for Q < m < k 

|^^+m^™(g-tr,(r-))| ^ q aS T ^ 0. 

Using the estimate (16) and the condition (13) we also have 

|^^+m^m^g-tr,(r))| < ^U(z)+o,va ^--tn(r) ^ q aS r ^ OO 

and therefore Claim 1 is proved. 

Claim 2: The following Mellin's inversion formula holds for c > 

Proof of Claim 2: If we show that Mt(z) is integrable, then by Mellin's inversion 
theorem (1) we get 

^-tvir) = J-[ Mtiz)r-'-dz 
27rj 
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and so we get the result from Claim 1. Using the estimate (16) we have 
\Mt^iz)\ < f |£l'^(e-*''(''))r"+'^-i|dr 



< kAt / r5f(2)-i+"e-*')('-)dr 

+ k{At)'' [ r5R(z)-l+afeg-t^(r)^^^ 

The integrals on the right converge absolutely because of the condition (13). We 
notice that for 3?(z) > —a the function M^{z) is bounded and holomorphic on the 
bounded strip. Therefore for fc > 2 

(18) \Mt(z)\ = (z) < ^^-^^ 

^ ' I n 71 r{z + k) ' ^ 1 + |3(z)|'=^ 

and hence A/((z) is integrable. 

Claim 3: Then K'^{t, x) has an integral representation for ((i+l)/2+/3 < c < d+P 

(19) (-1)' 1 / r(z)r(^)2-+^ 

Proof of Claim 3: Recall the identity (12), that is 

nOO 

if'3(t,a;) = (27r)-'^/2|^|-d/2+i / j^/2_^(|a;|r)r'i/2+/3g-«^W^^^ 

Jq 

By Claim 2 we obtain 

We now recall that for > —1/2 and r > 

(20) |J.(r)| <(r''Ar-i/2). 

We choose ((i+l)/2 + /3<c<(i + /3so that the integrand 

^lM_Mnz)J,/,-A\^\r)r~^'''/'^ 

is integrable from the estimates (18) and (20). So we can apply Fubini's theorem 
to obtain 

Since we have 

nOO 



a; 



d/2 p 1 2~-'+'^^^+^r( ~^+'^+'^ ) 



r(^) 

using the identity (2), Claim 3 follows. 
Claim 4: As \x\ — >■ oo, 

(21) j^^(t,a:)^ ^ ' 
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Proof of Claim 4: In the Proof of Claim 2, we noticed that the fmiction Mf{z) is 
holomorphic for 3fJ(z) > —a. Since r((— z + d + f3)/2) has no pole for < d + /3 
and l3 is not an even integer, the integrand 

V(z)r( -^+'^+/^ ^12-^+^ 

in the equation (19) is meromorphic in 5R(z) > —a and has a simple pole at z = 0. 
We calculate the residue at z = as follows: 

r(A;)r(^) r(fc)r(^) Jo 



r(^)2^ 



^ V 2 



^ ^ , (-1)^6-*"^") 



2 ^ 

by integrating by parts and using the Fundamental Theorem of Calculus. Now we 
set 

5= (aAl)/2 

and define the contour 7^ consists of a rectangle with vertices c — iR, c + iR, — c + 
iR, —c — iR with the counterclockwise orientation. Since the only pole of the inte- 
grand (22) inside the rectangle 7^ is at 2; = 0, Cauchy's Residue theorem gives 

1 f ^(z)^(^)2-^+^ (-i)^r(^ 



; a Mf(z)\x\^dz ^ —f^ — e~ 

By Stirling's formula, the two integrals over the horizontal segment on top and 
bottom tend to as i? tends to infinity. Thus, we can shift the vertical line left to 
obtain 

2'3r(^) e-*''(o) 



where 



7r'i/2r(^) \x\ 
{^lf2^ 1 f T{zn=^)2-^+\.^ 



R{t,x) ^ / ' ' ' ^ S Mt^{z)\x\^dz. 



Using Stirling's formula, we have the estimate 

r(z)r( ~"+/^^ )2-"+'3 
r(z + fc)r(^ 

and therefore \x\'^^^R{t,x) tends to as |a:;| tends to infinity. □ 



-Mnz)\x\^ < |cj(^)|-SR(.)+./2+/.-., 



Remark 6. The key point of the proof is that after inserting the Mellin inversion 
formula (17) to (12) and changing the order of integration, the remaining integrand 
(22) has meromorphic continuation. 

Remark 7. If 1 < a < 2, we choose c ~ (a + l)/2 so that the integrand (22) has 
two simple poles at 2: = and z = — 1. Then the residue at z = — 1 produces the 
second order approximation. Actually, one can go further, if the symbol 77 satisfies 
the regularity condition for all fc g N, then we can obtain an asymptotic expansion 
of the kernel. 
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Remark 8. The assumption (14) is analogous to that of standard symbol class in 
the theory of pseudo-differential operators. 

Remark 9. The meromorphic continuation of the Mellin transform using resolution 
of singularities goes back to Bernstein-Gelfand [BGG9]. 

We note that the main contribution to K^(t, x) comes from the integration near 
r = 0. We guess that the assumption (14) could be weakened away from the origin. 
This idea is made precise in the following proposition. 

Let ipir) be a smooth function defined in M, with the properties that ip{r) = 
for r < 1, and V'(^) = 1 for r > 2. Let us set for /3 > 

/oo 
Jd/2-i{\^\rmr)r''/'+^e-'^^-Ur. 

We shall show that this is an error term. 

Proposition 1. Suppose thatri(r) G C^((0, cx))) satisfies the integrahility condition 
(23) r r'-'^-^'^+^e-'^'^'-Ur <oo. 



Then for fixed < > 

\£^{t,x)\<\x\-^-''\ 
Proof. We begin by recalling that for v > —1/2 and r > 
(24) i?KJ,(r)]=r''J,_i(r), 
where D denotes the derivative with respect to r. Fix t > and denote 

g{r) =iP{ry/'^+'^e-*'''^''l 

Then using (24) we have 

J^-ii\x\r)g{r)dr 



oo 

{\x\rrj,M\x\r){\x\r)-''g{r)dr 

OO -I 

, -D[{\x\rrM\x\rmx\r)-''g{r)dr 


hm lj,(\x\r)g{r) - I H {\x\rr .M\x\r)D[i\x\r)"' g{r)]dr 

r-i-oo k ^ Jo 

r / ^ , s 1 



I Ji^(\x\r)r [vl — rD)g{r)dr 
k ' 

since lim^-i-oo Ji^i\x\'i')gif) — 0. If we denote 

L,=r-\iyI^rD), 

then we can write 

1 f°° 

J^-i{\x\r)g{r)dr ^ — / Jt,{\x\r)L^g{r)dr. 
m Jo 

Carrying out the repeated integration by parts gives 
1 f°° 

£'^(t,x) ^ —TTT / J^+N~i{\x\r)L^+N-iLu+N-2- ■ ■ Lug{r)dr. 
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\x 


pV+1/2 




1 


\x 


pV+1/2 




1 


\x 


A'+l/2 




1 


\x 


A'+l/2 




1 



Using the decay of Bessel functions (20) and the integrabihty condition (23), we 
have 

1 /'°° 

^^it^^)\^TZJMTU^ r-^'^\L,+N-iL,+N-2 - ■ ■ L,9{r)\dr 
/ r~^'^^r"-^\D''g{r)\dr 

■^1 n=0 

N „oo n 
< - V / r-l/2+n-JV^|^rx-fc^d/2+/3||^fcg-tr,(r)|^^ 

n=0'^l fc=0 



^ |a;|A+i/2 

and the result follows. □ 

We have the same main term under the weak assumptions compare to those of 
the previous theorem. This shows that the local behavior of the symbol 77 near 
r = is important. 

Corollary 5. Suppose that ri{r) is a real-valued continuous function on [0, 00) and 
that there exist < a < 2, k > (d + 3)/2 + /3, and M > such that 

(25) sup sup r-"+"'\D"'T]{r)\ < 00, 

l<m<k 0<r<l 

(26) sup sup iD'^rKOIr"^^ < 00, 

0<r?i<fe l<r<oo 



an 



d 



(27) liminf^^oo. 

1 — >oo log r 

If (5 ^ 0,2, A, ... , then for each t > 

K"(t,x) ^ : '^S--, — ,j,a 0,S \x\ — ?> 00 

Proof. By essentially the same way of in the proof of Theorem 2, we can prove the 
following Mellin inversion formula for c > 

271-1 J(c) ^[z + k) 

where 

M,'(=)=|^.«-(A)'e-*.<i„. 

The remaining proof is the same and so it is omitted. □ 

We consider the perturbation of the symbol corresponding to the symmetric 
a-stable process that has the form 77 (r) = r" + 771(7') under some conditions of 
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Theorem 3. Let < a < 2. Suppose that r]{r) = r" +rii{r) and that 771(7") is a real- 
valued continuous function on [0, 00) that satisfies for some 6 > a, k > {d+3) /2+/3, 
and M > such that 

(28) sup sup r-'+™|D™7/i(r)| < 00, 

l<m<k 0<r<l 



(29) sup sup |L»'"7;i(r)|r-^^ < cx), 

0<m</c l<r<oo 

and 

(30) liminf^^=(X). 

r— >oo log r 



If P = 0,2, A, ... , then for each t > 
K^{t,x) - - 



Proof. From the identity (12) we already know that 

(31) K^t,x) = (.^^./^i^,,^, [ J,,,_i(.)e-*(-/l-l)°e--''^(-/l-l)r'^/^+^dr. 
Since r(z) is the Mellin transform of e"'', we have for Cq > 

(32) f,-t{r/\x\r ^ 1 j r(-)(ii/"r/|x|)-^ds 



by Mellin's inversion formula (1). We also have for c > 

(33) e-..»<''N. = <^£fI<iL_Mf,,,„H)-.,,, 

where 

M^{z) = (ll) 

Putting (32) and (33) together with (31), we obtain 

K^it x) - i-^)''^^ r [ f r(2:)r(^) .^^ 

- a(27r)^/2+2|^|.+^ I y^^^^, J^^^ T{z + k) ^* 

X Jd/2-i{r)r-'^"^'^''^^'^dzdsdr. 
From the identity (2) we have 



00 



Jai.-.ir)r~^--^^l--^Ar = ^^^^ ^ \ 



18 



S. JO AND M. YANG 



We choose {d + 1) /2 + P < c + co < d + (3 so that we can apply Fubini's theorem 
and the identity (2) to obtain 

nOO 

X / Jd/2-i(r)r-'^+^+'*/2+'^drdzds 

"'0 



a(27r)'i/2+2|a;|'i+/3 4,,)4) r{z + k) 
^ dzds 



2 

(_;^)fc+i2<i/2+/? . . r(s/a)r(z)r( '"-"+'^+^ ) 

^ «(27r)'i/2+2|x|'^+/^ 4) r(z + fc)r(^±|^^ 



Now, we define the contour consists of a rectangle with vertices c — ii?, c - 
iR, —5/2 + iR, —(5/2 — iR with the counterclockwise orientation. The integrand 



r(z + fc)r(^±|^ 



Mf(z)(|x|/2)^ 



has only one simple pole at 2; = in the rectangle bounded by 7jj. By Cauchy's 
Residue Theorem, 

J^^r{z + k)r{^±^) r(fc)r(^) ' 

We have M/=(0) = (-l)'=r(A:)e-*''i(''l As i? ^ 00, the two integrals over the line 
parallel to x-axis go to zero by Stirling's formula. Thus, we can shift the line 
integral left to obtain 



2 



27rz7(,) r(z + fc)r(^+p^) r(^) 



where 



1 r r(z)r(^-^^^±^^) 

Rtx{s) = — ^-j—^M,^{z){\x\/2ydz. 
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Hence we have 

_2d/2+/3i -tr,i(0) r Y ( S / a)r ( 

^ ^ a(27r)<i/2+i|a;|''+^ r(^) ^' ^ 

+ a(27r)'^/2+i|x|'^+/3y(_3„) r(^) il^l/^)"!^ 

For > (rf+3)/2+;3, the last two integrals converge absolutely by Stirling's formula 
and they are bounded by 

|^|-d-/3-(5^^Afa)^ 

as \x\ — > cxD. Since { ^"^^ A |a) > a, the theorem is proved. □ 
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